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ABSTRACT 

We investigate the gravitational lensing properties of dark matter halos 
with Navarro, Frenk & White (NFW) density profiles and derive an analytic 
expression for the radial dependence of the shear, 7(0;), due to these objects. In 
addition, we derive an expression for the mean shear interior to a given radius, 
7(x), and use this to quantify systematic errors that will arise in weak lensing 
mass estimates of astronomical objects in the case that the mass estimate is 
based on an a priori assumption that the underlying potential is that of a 
singular isothermal sphere when, in fact, the potential is that of an NFW-type 
object. On mass scales between IO^^Mq < M200 ^ IO^^Mq, the assumption 
of an isothermal sphere potential results in an overestimate of the halo mass, 
and the amount by which the mass is overestimated increases linearly with the 
value of the NFW concentration parameter. Considerable overestimates of the 
mass (~ 60%) can occur for galaxy-sized halos, but for rich clusters the mass 
overestimate is small. The degree to which the mass is systematically in error 
is dependent upon the cosmology adopted, with a COBE-normalized standard 
CDM model yielding the largest systematic errors for a given true value of the 
halo mass. 

Subject headings: cosmology: theory — dark matter — gravitational lensing — 
galaxies: clusters: general 

1. Introduction 

Several recent numerical investigations (e.g., Navarro, Frenk & White 1997, 1996, 
1995) have indicated the existence of a universal density profile for dark matter halos that 
results from the generic dissipationless collapse of density fluctuations. Interior to the virial 
radius, the Navarro, Frenk & White (NFW) profile appears to be a very good description 
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of the radial mass distribution of simulated objects that span 9 orders of magnitude in 
mass (mass scales ranging from that of globular clusters to that of large galaxy clusters). 
The apparent generality of the NFW density profile has been confirmed independently by a 
number of studies (e.g., Bartelmann et al. 1998; Thomas ct al. 1998, Carlberg et al. 1997; 
Cole & Lacey 1997; Kravtsov, Klypin & Khokhlov 1997; Tormcn, Bouchet & White 1997); 
however, there are a few controversial claims that the NFW prescription may fail at very 
small radii (e.g., Ghigna et al. 1998; Moore et al. 1998). 

The NFW density profile is given by 

[r/rs) (1 + r/rs) 

where pc — is the critical density for closure of the universe at the redshift, z, of the 

halo, H{z) is Bubble's parameter at that same redshift, and G is Newton's constant. The 
scale radius = T2qq/c is a characteristic radius of the cluster, c is a dimensionless number 
known as the concentration parameter, and 

^ 200 

3 ln(l + c) - c/(l + c) ^ ' 

is a characteristic overdensity for the halo. The virial radius, r2oo; is defined as the radius 
inside which the mass density of the halo is equal to 200pc (see, e.g., Navarro, Frenk & 
White 1997). The mass of an NFW halo contained within a radius of r2oo is therefore 

N SOOtt s 8007rp(z) o 
M200 = M{r2oo) = -^Pcrloo = "§-^^200 (3) 

where p(z) is the mean mass density of the universe at redshift z and Q{z) is the density 
parameter at redshift z. 

Although it has not been proven categorically, it is certainly widely-thought that the 
masses of large galaxies, groups of galaxies, galaxy clusters, and supcrclusters arc dominated 
by some form of dissipationless dark matter. Therefore, it would not be unreasonable to 
expect that the spherically-averaged density profiles of these objects would be approximated 
fairly well by NFW profiles. Observationally, the total masses and mass-to-hght ratios 
of these objects are not constrained especially well at present; however, this situation is 
changing rapidly, due in large part to the fact that high-quality imaging of gravitational 
lens systems is yielding direct constraints on the nature of the mass distribution within the 
dark matter halos. 



Observations of gravitational lensing provide powerful constraints on both the total 
mass and the mass distribution within the lens itself, owing to the fact that one essentially 
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uses photons emitted by objects more distant than the lens to trace the underlying 
gravitational potential of the lens directly. In particular, large chistcrs of galaxies (which 
are both massive and centrally-condensed) are especially good gravitational lens candidates, 
and detections of the coherent pattern of weak lensing shear due to a number of clusters has 
led to interesting constraints on the masses of these objects (e.g., Tyson, Wenk & Valdes 
1990; Bonnet et al. 1994; Dahle, Maddox & Lilje 1994; Fahlman et al. 1994; Melher et al. 
1994; Small et al. 1994, 1995, 1997; Tyson & Fischer 1995; Small & Dickinson 1995; Kneib 
et al 1996; Seitz et al. 1996; Squires et al. 1996ab; Bower & Small 1997; Fischer et al. 1997; 
Fischer & Tyson 1997; Luppino & Kaiser 1997; Clowe et al. 1998; Hoekstra et al. 1998). 
Although more controversial than the results for lensing clusters, detections of systematic 
weak lensing of distant field galaxies by foreground field galaxies have been reported and 
these have been used to place constraints on the physical sizes and total masses of the dark 
matter halos of the lens galaxies (e.g., Brainerd, Blandford & Small 1996; Griffiths et al. 
1996; Ebbcls 1998; Hudson ct al. 1998; Natarajan et al. 1998). Additionally, a detection of 
the coherent weak lensing shear due to a supercluster has been reported recently (Kaiser et 
al. 1998). 

Because of the apparent direct applicability of the NFW density profile to the dominant 
mass component of all of these objects, and because of the potential of observations of 
gravitational lensing to provide strong, direct constraints on the amount and distribution 
of dark matter within them, we investigate the lensing characteristics of dark matter halos 
with generic NFW-type density profiles in this paper. In §2 we compute the convergence 
and the shear profiles of NFW halos. In §3 we compare the mean shear induced by 
NFW lenses to that of simpler singular isothermal sphere (SIS) lenses and consider the 
implications of our results for possible systematic errors in lens masses that are determined 
in observational investigations which invoke an a priori assumption of an isothermal lens 
potential. A discussion of the results is presented in §4. 



2. Convergence and Shecir of an NFW Object 

We perform all of our calculations below using the thin lens approximation, in which an 
object's lensing properties can be computed solely from a scaled, 2-dimensional Newtonian 
potential. The thin lens approximation is valid in the limit that the scale size of the lens 
is very much less than the path length traveled by the photons as they propagate from the 
source to the lens and from the lens to the observer. In this case the lensing properties of 
an object are completely described by two quantities, the convergence, /t, and the shear, 
7. The names of these quantities are indicative of their effects upon a lensed image; the 
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convergence describes the isotropic focusing of light rays while the shear describes the effect 
of tidal gravitational forces. Convergence acting alone leads to an isotropic magnification 
or demagnification while the shear induces distortions in the shapes of lensed images. 

If we define z to be the optic axis, then for a lens with a 3-dimensional potential 
^{DfiO, z) we can formulate a conveniently-scaled potential as projected on the sky: 

Here is a radius vector on the sky and D^, Dg^ and Dds are, respectively, the angular 
diameter distances between the observer and the lens, the observer and the source, 
and the lens and the source. Under the definition of ip{9) above, the convergence and 
the components of the shear tensor may be written as straightforward combinations of 

— * 

second-order derivatives of ^0 with respect to image plane coordinates 6 = (^i, ^2), 
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The magnitude of the shear is simply 7 = I7I = Y7i + 72 (e.g., Schneider, Ehlers & Falco 
1992). In the limit of weak gravitational lensing, the convergence and shear are formally 
small (i.e., k « 1, 7 << 1), the ellipticity induced in the image of an intrinsically circular 
source due to lensing is of order 7/2, and the position angle of the lensed image ellipse is of 
order the phase of 7 (e.g., Schramm & Kayser 1995; Seitz & Schneider 1997). 

The local value of the convergence may be expressed simply as the ratio of the local 
value of the surface mass density to the critical surface mass density: 

<S) = ^, (8) 

where 

(e.g., Schneider, Ehlers & Falco 1992) and c in the equation above is the velocity of light. 
The radial dependence of the surface mass density of a spherically symmetric lens such as 
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an NFW lens is obtained simply by integrating the 3-dimensional density profile along the 
line of sight, 

poo 

E{R) = 2 / p{R,z)dz, (10) 



where R = DayOf + 02 is the projected radius relative to the center of the lens. 

For convenience we will adopt a dimensionless radial distance, x = R/vg. Integrating 
equation (|l|) along the line of sight, the radial dependence of the surface mass density of an 
NFW lens can then be written as: 



Snfw(a;) = < 



2rs5cPc 
(x^-l) 



2rs5cPc 
3 



1 - 



VT- 



Farctanh 



l~x 
1+x 



2raScpc 



— / ^ arctan \ / — 



(x < 1) 
(x = l) 
(x > 1) 



(e.g., Bartelmann 1996). The radial dependence of the convergence due to an NFW lens is 
then simply Kntwix) = Snf„(x)/Sc. 

Since the NFW density profile is spherically symmetric, the radial dependence of the 
shear can be written as _ 

^nfwix) ^nfw(3^) 



7nfw(a;) 



(12) 



(e.g., Miralda-Escude 1991) where Snfw(a;) is the mean surface mass density interior to the 
dimensionless radius x. In terms of this radius, then, the mean surface mass density of an 
NFW halo is given by 



Snfw(3^) — / ^ Snfw(3^ )dx 

Jo 



^^sScPc [^arctanh^i^f + In (f ) 



4r,4p, [l + In (i 



^2 '^s^cPc 



arctan ^/T-rz + In ( | 



.ys^ Y i+x 

and the radial dependence of the shear is, therefore, 

'^^<(x) {x < 1) 



Sc 



ft + 4 In (1) 

9>{x) 



(x=l) 
(x > I) 



{x < 1) 

{x = 1) (13) 
(x > 1) 



(14) 
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where the functions (7<,>(a;) above depend upon only the dimensionless radius x and are 
exphcitly independent of the cosmology: 



Sarctanhj^-^ 4 /r 
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(15) 
(16) 



Equation (14) above can also be obtained straightforwardly from equations (7) through (11) 
of Bartelmann (1996). The radial dependence of the shear due to an NFW lens is shown in 
Fig. 1. 

The shear due to a given lens (e.g., a cluster of galaxies) is computed directly from 
the coherent distortion pattern that it induces in the images of distant source galaxies. 
In the realistic observational limit of weak shear and a finite number of lensed images, 
a measurement of the mean shear interior to a radius x centered on the center of mass 
of the lens (i.e., 7(0;)) is more easily determined than the differential radial dependence 
of the shear (i.e., 'j{x)). In the case of the NFW profile, the mean shear interior to a 
(dimensionless) radius x can be computed directly from equation (|14|) above: 



/O l^r 



X^ 



— g^{x')x'dx' + gy{x')x'dx' 



(17) 



A useful fiducial radius interior to which one might measure the mean shear is the virial 
radius, R = r2oo, or equivalently, interior to a; = (r2oo/''^s) = c, where c is the concentration 
parameter. For all masses of astrophysical interest c is greater than 1 and, therefore, the 
mean shear interior to the virial radius becomes 



7nfw 1^^200) 

which we rewrite as 



rsScPc 



g^{x)x'dx + / gy{x)x'dx 



7nfw('^2oo) = ^^(c) 



where 



— (^J g^{x')x'dx' + j gy{x')x'dx' 



is a function of the concentration parameter alone. 



(18) 

(19) 
(20) 
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3. Comparison to the Singular Isothermal Sphere 

Like the NFW mass profile, the singular isothermal sphere (SIS) mass profile is 
characterized by a single parameter (i.e., the velocity dispersion, a^). The mass of an SIS 
interior to a three dimensional radius r is: 

M{r) = (21) 

(e.g., Binney & Tremaine 1987) and the mean gravitational lensing shear interior to a radius 
R that is induced by an SIS lens is: 

7..(^) = ^^ (22) 

(e.g. Schneider, Ehlers & Falco 1992). 

Because of its simplicity, the SIS density profile is sometimes adopted in observational 
investigations in order to obtain an estimate of the mass of a lens without fully reconstructing 
its true underlying density profile (e.g., Tyson, Wenk & Valdes 1990; Bonnet et al. 1994; 
Small et al. 1994, 1997; Small & Dickinson 1995; Bower & Small 1997; Fischer & Tyson 
1997). By assuming that the underlying potential of the lens is well-approximated by an 
SIS, a measurement of the mean shear interior to a projected radius R leads directly to a 
measurement of the velocity dispersion of the lens (e.g., equation |2^), which in turn leads 
directly to an estimate of the mass of the lens (e.g., equation |2T| ). 

The NFW density profile, which is shallower than isothermal on small scales, and 
which turns over to isothermal on large scales has, however, been shown to be a far better 
approximation than the SIS to the spherically-averaged density profiles of halos formed via 
dissipationless collapse. Therefore, it is likely that lens mass estimates based on an a priori 
assumption of an isothermal potential will be systematically in error. In this section we 
compare the mean shear induced by NFW lenses to that induced by SIS lenses, under the 
constraint that the NFW and SIS lenses both have identical virial radii, r2oo; and, therefore, 
identical masses interior to r2oo- From this we will then investigate the possible systematic 
errors in lens mass estimates that would arise due to the assumption of an isothermal 
potential when, in fact, the lens is best represented by an NFW density profile. 

Let us consider two lenses which have identical masses, M200, interior to the virial 
radius. One of the lenses has an NFW density profile with a concentration parameter of 
c and the other is a singular isothermal sphere with velocity dispersion a^. If these two 
objects have identical redshifts, z^, and act as lenses for populations of source galaxies 
which have identical redshifts, Zg, then from equations (|I9|) and ( ^2]) above, the ratio of the 
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mean shears induced by these two lenses interior to r2oo is given by: 

7nfw ('"aoo) _ rsPcr200 , n 

— r- — LrJ- [C . 

7sis (^"200) 

Using equations (§) and (pT]) above and recalhng that the concentration parameter is 
c = r20Q/rs, it is straightforward to show that equation ( ^31) reduces to 



which is a function solely of the concentration parameter of the NFW lens and is 
explicitly independent of the redshift of the sources, z^. Because of the dependence of the 
concentration parameter on both the redshift of the lens and the cosmology through 'p{zd), 
equation ( plD is not explicitly independent of either the cosmology or the lens redshift, 
Zd- However, for lenses of a given mass, its dependence on both z^ and the cosmology is 
relatively weak. 

Shown in Figs. 2 and 3 are the ratio of the mean shears interior to r2oo for NFW and 
SIS lenses with virial masses in the range of IO^'^^Mq < M200 ^ IO^'^^Mq. Fig. 2 shows 
the results for lenses located at Zd = 0.1 and Fig. 3 shows the results for lenses located 
at Zd = 0.5. The four panels in the figures show the effects of varying the cosmology, and 
plotted along the top axes of all of the panels is the NFW concentration parameter which 
corresponds to the lens mass plotted on the lower axes. 

Two of the cosmologies illustrated in Figs. 2 and 3 are standard cold dark matter 
(CDM) cosmologies, which differ from one another only in the choice of the normalization 
of the power spectrum (SCDM-I is a cluster abundance normalization while SCDM-II 
is COBE- normalized). The other two cosmologies are an open CDM model with zero 
cosmological constant (OCDM) and a spatially flat, low matter density CDM model with a 
large cosmological constant (ACDM). The parameters adopted for each of the models are 
summarized in Table 1 where Aq = X/SH^, Hq = lOO/i km/s/Mpc, n is the index of the 
primordial power spectrum of density fluctuations and 



7nfw (^200) ^ 3 J^(C) 

7sis (^^200) 4007r c 



(24) 




(25) 



Table 1: Cosmological Model Parameters 
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^0 


Ao 


h 




n 


SCDM-I 


1.0 


0.0 


0.50 


0.63 


1.0 


SCDM-II 


1.0 


0.0 


0.50 


1.20 


1.0 


OCDM 


0.25 


0.0 


0.70 


0.85 


1.0 


ACDM 


0.25 


0.75 


0.75 


1.30 


1.0 



The FORTRAN program charden.f, written and generously provided by Julio Navarro, 
was used to calculate the values of the concentration parameters for the NFW lenses in the 
above cosmologies. For each of the cosmologies, c was determined for halos with masses in 
the range of IO^^Mq < M200 < IO^^Mq at redshifts of Zd = 0.1 and za = 0.5. These values 
of c were then used in conjunction with equation (|2^ ) to compute the ratio of the NFW to 
SIS mean shear interior to the virial radius. 



For a given cosmology, it is clear by comparing Fig. 2 with Fig. 3 that equation (p4D is 
only weakly dependent on the lens redshift, Zd- The largest difference between the various 
panels in Figs. 2 and 3 which correspond to identical cosmologies occurs for SCDM lenses 
with masses ~ IO'^^Mq, and in this case the difference between Zd = 0.1 and Zd = 0.5 is only 
~ 10%. Similarly, by comparing the results plotted in all of the individual panels of Fig. 2 
and Fig. 3 at fixed Zd, it is clear that equation ( p4|) is not tremendously sensitive to the 
cosmology. In particular, the ACDM, OCM, and SCDM-I models all yield functions with 
nearly identical amplitudes for a given value of Zd- The SCDM-II model yields a function 
which is somewhat higher than the other three models, exceeding the others by ~ 25% for 
halos with masses ~ IO^^Mq and by ~ 20% for halos with masses ~ IO^^^Mq. 

Over the majority of the mass range investigated, the NFW lenses give rise to a mean 
shear interior to r2oo which is systematically larger than that of the SIS lenses. As a result, 
if one were to measure the mean shear interior to a radius of r2oo of an NFW halo, yet 
assume it to be an isothermal sphere, the resulting estimate of the virial mass of the lens 



(M200) would be systematically high. From equations ([21| ) and p^ ) above, it follows that 
the mass of an SIS lens interior to r2oo is simply: 

M200 = 2Scr^oo7(r2oo) (26) 

so that the mass inferred for the lens scales linearly with the mean shear. Therefore, the 
systematic error in the true virial mass of the lens is simply the ratio of the mean shear due 
to an NFW lens to that of an SIS lens with an identical amount of mass contained inside 
r2oo (i-e.. Figs. 2 and 3). 

Shown in Fig. 4 is the ratio of the mean shear (interior to r2oo) of an NFW lens and 
an SIS lens, plotted as a function of the NFW concentration parameter. (As in Figs. 2 and 
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3, both lenses have identical masses interior to r2oo)- From this figure, then, if one were to 
measure a mean shear for a given NFW lens, yet model the lens as an isothermal sphere, 
the degree of systematic error in an estimate of the virial mass would clearly be a function 
of the concentration parameter of the lens. For a given halo mass, the concentration 
parameter is a function of the cosmology (e.g., the top axes of Figs. 2 and 3); however, it is 
always the case that for a given cosmology, the larger the value of c, the lower is the value 
of M200- The general conclusions that can be drawn from Fig. 4 are: [1] the lower the mass 
of an NFW halo, the larger the systematic error in the mass estimate if the lens is assumed 
to be an isothermal sphere and [2] for a halo of a given mass, the largest systematic error 
in the mass estimate occurs in a COBE-normahzed cosmology (i.e., SCDM-II). With the 
exception of SCDM-II for which the error is somewhat larger, the systematic error in an 
estimate of M200 for rich clusters (M200 ~ IO^^Mq) is negligible (< 10%). The systematic 
error in an estimate of M200 for galaxy-mass objects (M200 ~ IO^^Mq) is, on the other 
hand, considerable (of order 55% to 65% for the SCDM-II model and of order 30% to 40% 
for the other models). 

4. Discussion 

It is generally thought that the masses of large galaxies, galaxy groups, galaxy clusters, 
and supcrclustcrs arc dominated by some form of dissipationless dark matter and, thus, it 
is not unreasonable to expect that their underlying mass density profiles will be represented 
reasonably well by NFW profiles. In addition, since their total masses and mass-to-light 
ratios are not strongly constrained at present, a significant effort is currently being devoted 
to the use of observations of gravitational lensing by these objects to quantify the amount 
and distribution of dark matter within them. We have, therefore, investigated the properties 
of NFW lenses in this paper and we have presented analytic expressions for the radial 
dependence of the convergence, k.{x), and shear, j{x), due to dark matter halos which have 
NFW density profiles. We have also presented an expression for the mean shear interior to 
a given radius, 7(x), due to NFW lenses and we have compared the mean shear interior to 
the virial radius of an NFW lens to that yielded by a singular isothermal sphere lens with 
an identical virial mass. 

It is not uncommon for the mass of a gravitational lens to be estimated under an 
assumption that the lens may be approximated by a singular isothermal sphere. However, it 
has been clearly demonstrated that the NFW density profile is a far better approximation 
to the density profile of objects formed by generic dissipationless collapse than is the 
isothermal sphere. We have computed the systematic error that would be encountered 
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in an estimate of the mass of an NFW lens, where the lens is assumed a priori to be an 
isothermal sphere. Over mass scales of IO^^Mq < M200 ^ IO^^Mq, the mass of the NFW 
lens is systematically overestimated when it is assumed that, for a given measured value of 
7(^200)) the lens can be approximated by an isothermal sphere. 

The size of the systematic error in the lens mass due to the isothermal sphere 
assumption is a function of the NFW concentration parameter of the lens, with the largest 
error occurring for halos with the largest values of c and, hence, with the smallest masses. 
The systematic error in the mass is not dramatic (i.e., not even as much as a factor of ~ 2), 
but this is unsurprizing since the shape of the NFW density profile in the outer regions of 
the halo is fairly close to an isothermal profile. 

In the case of halos with masses comparable to that of rich clusters, M200 ~ IO^'^^Mq, 
the systematic error in the mass due to the assumption of an isothermal potential is small. 
Therefore, the masses of lensing clusters that are estimated under the assumption of an 
isothermal potential (and in the limit that the shear is detected out to a radius that is large 
enough to be comparable to r2oo) should not have large systematic errors if, indeed, their 
density profiles are fitted well by NFW profiles. 

However, recent observations of lensing of distant field galaxies by nearby field 
galaxies (and, additionally, by the individual galaxies within clusters, e.g., Natarajan et 
al. 1998) have inspired a number of investigations through which the mass and extent of 
the dark matter halos of the lens galaxies might be constrained. The technique, known 
as galaxy-galaxy lensing, seems very promising at the moment, and in the near future a 
considerable amount of effort will be devoted to the use of observations of galaxy-galaxy 
lensing to constrain the nature of the dark matter halos of galaxies. The results of 
our investigation of systematic errors in the mass estimated for NFW lenses under the 
assumption of an isothermal potential indicate that these errors can be significant for 
galaxy-mass lenses (~ 60% in the case of a COBE- normalized CDM universe). Therefore, 
in the upcoming studies of galaxy-galaxy lensing, should an observational constraint on the 
masses of galaxy halos be based upon the assumption of an isothermal potential, it will 
be important to keep such systematic errors in mind when judging the strength of such a 
constraint. 
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Fig. 1. — Solid line: the dependence of 7nfw on the dimensionless radius x = R/vg, i.e., 
equations ([15D and (pIBI). Also shown for comparison (dashed line) is the radial dependence 
of the shear due to a singular isothermal lens, 7sis(a^) (x x~^. The normalization is arbitrary. 
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Fig. 2.— Ratio of the mean shear interior to a radius of r2oo 

for NFW and SIS lenses. 

Both lenses are constrained to have identical masses interior to r2oo (see text). The panels 
correspond to four different CDM models, and in this figure all of the lenses were placed at 
redshift of = 0.1. 
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Fig. 3. — Same as Fig. 2, but the lenses were placed at a redshift of — 0.5 
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Fig. 4. — Ratio of the mean shear interior to a radius of r2oo for NFW and SIS lenses as a 
function of the halo concentration parameter. The point types refer to four different CDM 
models. The left panel shows the result obtained when all lenses are placed at a redshift of 
Zd — 0.1; the right panel shows the result obtained when all lenses are placed at a redshift 
of Zd — 0.5. 



